Statistical mechanics of dissipative systems 
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We propose a generalization of classical statistical mechanics which describes the behavior of dissipative 
systems placed in contact with a heat bath. In contrast to conventional statistical mechanics, which assigns 
probabilities to the states of the system, the generalized theory assigns probabilities to the trajectories of the 
system. The conditional probability of pairs of states at two different times is given by a path integral. We 
present two simple analytically-tractable examples which illustrate the predicted effect of temperature on the 
mean trajectories, hysteresis and drift of the system. 
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The fundamental question addressed in the present 
communication is: How does a dissipative system behave 
when it is placed in contact with a heat bath ? It should 
be carefully noted that the systems under consideration 
are dissipative ab initio, e. g., as a result of internal fric- 
tion, viscosity, or some other dissipative mechanism. In 
particular, the systems are irreversible, path dependent, 
and exhibit hysteresis. Therefore, the objective of the 
theory is not to understand how dissipation arises from 
the coarse-graining of a conservative system or in con- 
nection with first-order phase transitions (see [[j], ^, [}| Q| 
for notable examples of that line of inquiry). Instead, 
the focus is in understanding how a heat bath influences 
the behavior of a dissipative system; what are the statis- 
tical properties of its trajectories; and how the effective 
behavior depends on temperature. 

For definiteness, we consider systems whose state is 
defined by an TV-dimensional array q of generalized co- 
ordinates. The energetics of the system is described by 
an energy function E(q,t). The explicit time depen- 
dence of E may arise, e. g., as a result of the application 
to the system of a time-dependent external field. In ad- 
dition, the system is assumed to possess viscosity, and, 
thus, the equilibrium equations are of the form: 



dE 
dq 



q,t) + f vis (q)=Q 



(1) 



where / Vls are the viscous forces. These equilibrium 
equations define a set of ordinary differential equations 
which, given appropriate conditions at, e. g., t — 0, can 
be solved for the trajectory q(t), t > 0. 

If the system is conservative, i. e., if / VIS = 0, the 
instantaneous state q(t) of the system at time t follows 
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directly from energy minimization, i. e., from the prob- 
lem: 



E(q(t),t) 



min E(q't) 



(2) 



In order to extend this variatonal framework to dissipa- 
tive systems, we resort to time discretization, leading to 
a sequence of minimization problems of the form (^|) 
[|5|, ^, [7p. Thus, we consider a time-discretized incre- 
mental process consisting of a sequence of states q ra _i 
at times to = 0, ... , t n -\, t n — t n -\ + At, .... We 
additionally assume that the viscous forces / Vls derive 
from a kinetic potential <fi(q) through the relation: 



and introduce the incremental work function: 

W(q n ,q n -i) = E(q n ) - £7(q„_i) 



mm 

{paths} 



(q(t)) dt 



(3) 



(4) 



where the minimum is taken over all paths q(t) such that 
<?(*n-i) = Qn-i an d <l{tn) = Qn- The fundamental 
property of the incremental work function W(q yl , q n -i) 
is that it acts as a potential for the forces f n at time t n , 
i. e., 



. _dW. 

Tn — "o — 9m Qn-lJ 
aq n 



(5) 



In order to verify this property, we may consider a small 
perturbation q n — > q n + Sq n , leading to a correspond- 
ing perturbation of the path q(t) — > q(t) + 6q(t), with 
8q(t n -i) = 6q n -i = 0. Then, in follows that 



5W(q n ,q n -i) = —(q n )-5q n + f™ 



f vis {q(t))-6q(t)dt 



5q n 



(6) 
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But the integral on the right hand side of this equation is 
to be evaluated along the minimizing path, and hence it 
vanishes identically. This gives the identity 



/™ • Sq n (7) 



QE 

SW(q n ,q n -i) = <{ -g-(q n ) 



Since the variation Sq n is arbitrary, it follows that 

dW dE 

^-(g„, g„-i) = g-(q n ) + /r = /„ (8) 

as stated. From this property it follows that the equilib- 
rium equation /„ = is the Euler-Lagrange equation 
corresponding to the minimum principle: 



minVF(q„,q„_i) 

<3n 



(9) 



This behavior is indistinguishable from that a conser- 
vative system with 'energy' W(q n) q n -i)- However, it 
should be carefully noted that W is determined by both 
the energetics and the kinetics of the system. Conse- 
quently, W depends on the initial conditions for the 
time step and, therefore, varies from step to step, which 
in turn allows for path dependency and hysteresis, as re- 
quired. A simple situation arises when the kinetic poten- 
tial 4>(q) is convex and coercive, i. e., it grows as |<ji| p , 
for some p E (1, oo), for large \q\. Then, the minimizing 
path is unique and consists of a straight segment joining 
q n -i and q n . Under these conditions, the corresponding 
incremental work function is 



W(q n ,q„-i) = E(q n ) - E(q n -i) + At<f> 



m 



(10) 



where we write Aq = q n — q n -\- 

Now, imagine placing the system just described in 
contact with a heat bath at absolute temperature T. The 
objective is to predict the ensuing behavior of the system. 
For conservative systems, such behavior is predicted by 
standard Gibbsian statistical mechanics, to wit: the prob- 
ability of finding a state q at time t obeys Boltzmann's 
distribution, i. e., is proportional to exp{— j3E(q, t)}, 
where (3 = 1/fceT and ks is Boltzmann's constant. In 
particular, the states of the system at two different times 
are uncorrected. We proceed to show that the analogy 
between the minimum principles (^) and (Jfy provides a 
basis for generalizing Gibbs's prescription to dissipative 
systems. However, in contrast to conventional statistical 
mechanics, which assigns probabilities to the states of 
the system, the extended theory assigns probabilities to 
the trajectories of the system. 

We begin by assuming that the incremental processes 
are Markovian, i. e., q n is correlated to q n -i but not to 
earlier states. Letp(q„, q n -i) be the joint probability of 



q n and q n -\, and let 

P(q n -l) = J p(q n ,q n -i)d N q n (11a) 

P(<ln) = J p(q n ,q n -i)d N q n -i (lib) 

be the corresponding probabilities of q n and q n _i. Then 
the conditional probability of q„ given q n -i is 



p{q n \qn-i) 



p(g«,g»-i) 
p(qn-i) 



(12) 



This probability may also be interpreted as the transition 
probability from state q n -i to the new state q„. We pos- 
tulate that the transition probability p{q n \q n -i) is Gibb- 
sian, i. e., 



p{q n \q n -i) 



i 



Z{q n -x) 



-/3W(q„,g„_i) 



(13) 



where 



Zfa-l) = J e -W*»*»-0 d N q n (14) 

is an incremental partition function. It therefore follows 
that: 



P(Qn 



1 



Z(q n -i) 



-/3W(q„,q„_i) 



p(q n -i) d N q n -i 



(15) 



Iterating this relation we obtain 

p(q n 




II z{ 1 P(lo) d N q n -i ■ ■ ■ d N q 

(16) 



In the limit of At — * and n — * oo at fixed t = nAt, 
( p^ ) defines a path integral of the form (see, e. g., [||]): 



p(qt\qo) 



-/3W[q] D 



(17) 



for some system-dependent work functional W[q]. This 
path integral gives the conditional probability p(qt\qo) 
of finding the system in state q t at time t > given that 
the system is in state qo at time t = 0. The integrand of 
( |l7| ) may be interpreted as relating to the probability of 
individual trajectories q(r). More precisely, given two 
trajectories qx{r) and q% (r), their relative probabilities 
are: 



pm\ 



-pw[ qi 



p[q 2 ] 



(18) 
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FIG. 1: Evolution of the probability density for a spring and 
dashpot system subjected to time-proportional loading. 



It is noteworthy that, in contrast to the conservative case, 
the presence of kinetics renders states of the system at 
different types correlated. However, this correlation may 
be expected to decay with elapsed time, conferring the 
system a. fading memory property. 

In the remainder of this letter we present two simple 
analytically-tractable examples which provide a first il- 
lustration of the theory, namely: a linear spring-dashpot 
system; and dry friction. In addition to their value as 
illustrative examples, these simple models should also 
useful as a basis for constructing mean-field approxima- 
tions to more complex systems. The treatment of general 
systems requires the use of numerical schemes such as 
the Path-Integral Monte Carlo (PIMC) method [§]. 

The linear spring-dashpot system is characterized by 
an energy function and kinetic potential of the form: 



E(q,t) = -q 2 -f(t)q 



(19a) 
(19b) 



where q £ R, k is the spring constant, r\ is the dashpot 
viscosity, and / (t) is a time-dependent applied force. For 
this system, the equilibrium equations (|l|) reduce to: 



Kq — f + rjq = 



(20) 



The corresponding incremental work function for this 
system is: 



W(q n ,q n -x) = <?n_i) 

— f n q n + Jn-Xln-l + (in — Qn-l) 



(21) 



pie Gaussian integral. Inserting the result in ( |l6| ) yields 

-n/2 

PiQn) 



(3{k + 77/ At) 



cxp < -/3 ^2 

{ i=0 

p(qo)dq n -i ...dqo 



(At (/j+i - Kq i+1 ) - 7](q i+1 - q t )Y 



2At(r) + nAt) 



(22) 



In the limit of At — > and n — > 00, with t = nAt fixed, 
the sum in the exponential of this formula converges to 
the Riemann integral 

W[q}= [ ±- (nqir) - f(r) + V q(r)) 2 dt (23) 



In the same limit, ( J22| ) may be written as the path integral 

\ -n/2 



Dq= lim 



2tt 



(3(K + r]/(t/n)) 



dq n -\ ■■■dq 
(24) 



which is of the anticipated form (|17|). The structure of 
( p^ is revealing. Thus, it is observed that the work func- 
tional W[q] integrates in time the square of the equilib- 
rium equation (EQ), reduced to units of power by means 
of the factor 1 pq7 Evidently, the path which contributes 
the most to the path integral is the critical path, i. e., 
the path q c which minimizes the work functional (|23|). 
Indeed, in the limit of zero viscosity or zero tempera- 
ture, the only path which contributes to the path integral 
is the critical path, which in that limit coincides with the 
deterministic trajectory, i. e., with the solution of (|20|). 
For finite viscosity and finite temperature, however, all 
paths contribute to the path integral to varying degrees, 
the contributions becoming increasingly weaker as the 
trajectories depart from the critical path. The stationar- 
ity of W[q] yields the Euler-Lagrange equations for the 
critical path 



w 2 q 



f_ 

V 



(25) 



and the evaluation of the partition function entails a sim- 



where u> = k/t) and the solution q c is subject to the 
boundary conditions q c {0) = qo, Qc(t) = qt- It is in- 
teresting to note that this equation is of second order in 
time, whereas the original equilibrium equation (|o|) is 
of first order. This increase in order makes it possible to 
enforce boundary conditions at times and t, in contrast 
to the original equation (^0|) which only allows for initial 
conditions. Furthermore, we note that (^5|) is obtained 
by squaring (^o|), and thus the solutions of the latter are 
subsumed within the former. However, by virtue of its 
higher order the Euler-Lagrange equation (^5j) has solu- 
tions which do not satisfy (20h. It thus follows that the 
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FIG. 2: Trajectories of a system sliding under the action of 
an applied cyclic load f(t) — A sin(wi) against a frictional 
resistance s. 



critical path does not coincide with the deterministic tra- 
jectory in general. Conveniently, the path integral ( |l7| ) 
can be evaluated exactly in closed form (e. g., and 
0]). The result is: 



2-7T sinh(o)i) 



-0{V[Qc]) 



(26) 



where the critical path q c (t) is the solution of (^). Fig. [j] 
shows the evolution of the probability density for a lin- 
ear loading history f(t) = At in terms of the nor- 
malized variables q — q/(Arj/K 2 ), t — t/(r]/K), and 
f = T/(A 2 ri 2 /k B K 3 ). The mean value (q(t)) traces the 
critical path q c (t), which, as noted earlier, differs from 
the classical path at finite viscosity and temperature. An 
additional effect of viscosity, which is clearly evident in 
these figures, is to cause the width of the probability den- 
sity to broaden in time, with the attendant increase in the 
uncertainty of the state of the system. 

The work functional ( |23| ) appearing in the path integral 
representation ( p^| ) is characteristic of systems whose ki- 
netic potential <fr(q) exhibits quadratic behavior near the 
origin. However, there are cases of interest which do not 
fall into that category. A case in point is provided by 
dry friction, which is characterized by a kinetic poten- 
tial which has a vertex at the origin. Consider, by way 
of example, a one-dimensional system sliding against a 
frictional resistance s > under the action of an applied 
force f(t) such that \f(t)\ < s. In this particular case 



E(q,t) = -f(t)q 
4>{q) = s\q\ 



(27a) 
(27b) 



and the incremental work function (|j) specializes to 

W(q n , q n -i) = -f n q n + fn-lQn-l + s \q n - ?n-l| 

(28) 



As expected from the rate-independent nature of dry fric- 
tion, W is independent of At. The corresponding parti- 
tion function ([m]), transition probability (pj[), and proba- 
bility density ( jig ) take the form: 



Z(q n -l) 



2s 



o/3(/n-/n-l) 3 n-l 



P(s 2 ~ fl) 



(29) 



(30) 



and 

p(q n ) = 



2s 



n 



cxp (s - Qi\ ~ - qi)) | 

p(q ) dq n -x . . . dq 

(31) 

respectively. An important feature of the transition prob- 
ability ( p0| ) is that it depends solely on the difference 
qn — QVi-i- Therefore, the chain qo, q%, . . . , q n defines a 
random walk, and the probability p(qt\ qo) corresponding 
to the limit of At — > and n — > oo at constant t = nAt 
follows by an application of the central limit theorem. 
The result is: 



p(qt\qo) 



l 



'(q t -{q)(t)) 2 /2a 2 (t) 



2ira(t) 



(32) 



where 



n-l 



(q){t)=qo+ lim - V 



2fi 



™ n U ^ ~ fi) (33) 
1 f 2/(r) 

is the mean path traced by the system, and 

2(s 2 + f?) 



a 2 (t)= lim ±Y1 

71. — >r>Ci 71 < " 



2( S 2 + f(r)) 



(34) 



measures the deviation from the mean path. It is in- 
teresting to note that, in the limit of zero temperature, 
(q)(t) = qo for as long as \f(t)\ < s, with the system 
sliding off to ±oo instantaneously as soon as / reaches 
±s. By way of sharp contrast, at finite temperature the 
system undergoes sliding even when the applied force 
fit) remains strictly below s in magnitude at all times. 
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In addition, the standard deviation a from the mean path 
is predicted to be proportional to temperature. By way 
of illustration of this behavior, the trajectories of the sys- 
tem corresponding to applied cyclic loads of the form 
f(t) — A sm(ujt) are shown in Fig. ^ in terms of the 
normalized variables: q = qj (ksT/s), f = f/s, and 
i = uit. The ability of the theory to allow for hysteresis 
and to predict its temperature dependence is particularly 
noteworthy. As expected, the response of the system is 
predicted to soften with increasing temperature, in keep- 



ing with numerical simulations of magnetic systems ni 
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